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1 .0 SUMMARY 


This report describes the modifications and additions made 
to the Flexible Spacecraft Dynamics (PSD) Program under contract 
NAS5-25678. The principal addition to the program was the cap- 
ability to simulate the Dynamics Explorer ~B Control System. The 
formulation for this addition is given in Section 2.0 of this re- 
port. The details of the modifications made to the PSD Program 
are presented in Section 3. Modifications to existing subroutines 
are briefly described and a detailed description of new subrou- 
tines is given. In addition, the program variables in new labelled 
COMMON blocks are described in detail. Section 4 gives a descrip- 
tion of new and modified input and output for the PSD Program. 


2.0 mathe mat ical fo r mulation f o r control system simulation 

A block diagram of the system to be simulated is given in 
Figure 2-1. This block diagram is represented as a set of first 
order ordinary differential equations, which are integrated in 
parallel with the equations of motion for the rest of the space- 
craft, using the same time step and integration algorithm. (The 
subroutine ADMIMP) . 

For the most part, the control components are linear dyn- 
amic systems. For such components, the stated transfer func- 
tions have been converted to state variable equations using stan' 
dard techniques. This transformation, however, is not always 
unique. Hence, it is necessary to state the exact form utilized 
in each case. 

In the equations following, the subscripts 1, 2 etc. are 
used primarily for convenience. However, the ordering of var- 
iables is the same an in program code. Hence a fourth order 
model with state variables x^...x^ may appear in program code as 
xg...Xg. The actual subscripts used in program code are given 
in the section on program inputs. Also, in this section, the 
symbols u and y denote (respectively) the input to and output 
from the given transfer function. In the system simulation, 
the blocks are coupled together. 
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Nutation Damper Phase Shift Circuit 
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y (the minus sign includes the sign inversion corres- 

ponding to path A in Figure 2-1) 

The inversion of sign at various roll angles and the switch be- 
tween primary and secondary damper may be simulated by changing 
the input parameters. 


Tachometer 

The transfer function 
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is represented as 
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For the reference control system, the tachometer significantly 
faster than the other dynamic elements. This causes the tach- 
ometer equation to dominate the time step control in numerical 
integration, which having little effect on system performance. 
Replacing the above transfer function with the static operator' 
y = K f u permits a significant reduction in program execution 
time. This altc .,tivc model is optionally available in the 
modified program. 
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Pitch Compensation Amplifier 


In unsaturated operation, the transfer function 

K e (i + s r t ) 

(i+ sTl) 

is represented as 

hh x < ' 

y= x t + Ko. •==• 1/ 

'z. 

Saturation occurs if / V/ ^ V^ )iy j 

If this occurs, Y is replaced by sign(Y), 


3.0 PR OGRAM A DD ITIONS AND MODIFICA T IONS 


The additions and modifications to the PSD Program accom- 
plished under this contract are described below . 

3.1 Summary of Modifications by Task 
Task 1 DYNAMICS EXPLORER-B CONTROL SYSTEM 

A simulation capability for the Dynamics Exploj B Con- 
trol System has been added to the PSD Program. 

1 . 1 Sensor Module 

The sensor module converts the state vector from the PSD 
simulation to the appropriate attitude angles (yaw, roll, and 
pitch) for sensor output. The pitch angle is compared to 
command pitch attitude and the resulting angular roll and pitch 
error are passed through either second order or fourth order 
lag circuits to the control system module. Provision has been 
made for adding bias errors and noise to the attitude angles ob- 
tained from the F8D state vector. 

1 • 2 Control System Module 

The control system module provides the driving torque to 
a variable speed momentum wheel. The driving torque is deter- 
mined from a simulation of the momentum wheel motor. The motor 
input is the sum of a constant bias speed voltage and a voltage 
derived from pitch and roll attitude error signals. The roll 
error signal is used for nutation damping. The pitch error sig- 
nal is used to maintain a constant pitch reference direction. 


The roll error signal is clipped and passed through a phase 
shifting circuit and the output is added to the pitch error 
signal* The resultant sum is passed through a compensation 
amplifier with voltage limiting. 

The print and plot subroutines of the PSD Program have 
been modified to permit printing and plotting of control sys- 
tem variables. 

The control sy s' -i transfer functions are integrated as 
added state variable 1 the PSD integrator ADMIMP. 

The DE-B control system simulation has been checked out 
with rigid body examples supplied by GSFC. t 

Task 2 EXTERN AT. MOMENTS 

Provision has been made to print and plot the sum of all 
external moments acting on the system. This sum includes gra- 
vity gradient, solar pressure, aerodynamic pressure, magnetic, 
and all thruster related control torques. 

Task 3 NOISE GENERATOR 

A Gaussian noise generator has been inclxaded so that the 
effect of noise on control system performance can be evaluated. 
This generator is based on subroutine GAUSS taken from the IBM 
System/360 scientific subroutine package. 

Task 4 TIIRUSTER OPTION 

The thruster option has been modified so that two thrust 
pulses can be applied to the spacecraft during a rotation period 


or time interval. The parameters defining the pulse character- 
istics are independent. 

Task 5 PAST FOURIER TRANSFORM OPTION 

The Fast Fourier Transform Option has teen modified so 
that up to four frequencies and amplitudes can be printed. 

Task 6 PROGRAM DISCREPANCIES 
6 . 1 Spin Axi s Moment Option 

The spin axis moment subroutines have been corrected. 

6 • 2 Deployment Acceleration 

The deployment acceleration has been reset to zero at the 
end of each deployment phase. 

6 . 3 Exces s ive X 0 Time 

The overlay structure has been adjxisted so that overlay 


10 time will be acceptable, 


3 .2 DKS CiKrVTION OF nUBROirJMN E MODI FIC WPIONS 

The subroutine modifications are briefly described below. 
ADMIMP 

Eliminate initialisation logic for integrator messages (moved 
to subroutine B0UNDS) » E3 iminntc common block JNTEGP. Call 
subroutine RWHOUT after every Integration step . 

MPNps 

Eliminate arguments in calling sequence. Establish default 
(large) values for bounds not specified via input. Initialise 
integrator message variables. 

CSD 

bogie changes so that the number of frequencies is specified by 
the input v.'ord I CSD. 

E CII0 A 

Echo printout was added for control system parameters. The out- 

» 

put for the thruster option was changed to reflect two pulses 
per spin period. 

GPRTNT 

Call to DEB0UT was inserted to provide for printed output of con 
trol system behavior. 

MAIN 

Modify calls to subroutine NUM, BOUNDS to eliminate arguement. 
Call subroutine CSIC, VDIC, Print program version number after 
every call to READ IN . Calls subroutine PDTAPR instead of DTAPRE 


Arrays DKR, dep now in labelled common AbfVTAT. References 
common block XCNTRL. Eliminate common blocks CNI50DY, VECTRS. 


NUM 

Eliminate arguments in call sequence. Calls NUMCSE if control 
simulation is specified. 

Modification for counting of pulses under conditions of two 
pulses per spin period. 

PULS ER 

Required changes to allow two pulses per spin period. 

PULS UH 

Required changes to allow two pulses per spin period. 

RD47R 

Setup calls added for control system simulation input. 

Setup calls and c0MMON/ f i’nRUST/changed for two pulses per spin 
period. 

sir 

Rex.' .3 the direction cosine matrix normalisation. Provide tran 
fer of rigid body angular rates from DEPEND array to OMEG vector 
(in common bl ock RPOOL1) . 

T0 TXMP 

Required changes to allow two pulses per spin period. 

T TRUS T 


Required changes to allow two pulses per spin period. 


wmyas 

Added control equations for closed loop control of momentum 
wheel speed. 

WRTPLT 

Plotting f ’ability added for body momentum vector components, 
external moments, and control system state vector. 
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! * 3. 3 DE SCRIPTION OE NEW SUBROUTINES 

The new subroutines written for the DE-B control system 
simulation are described in the folliwing section. 
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SUBR OUTIN E NAME : 


CSIC 


LANGUA GE : FORTRAN IV 


CALLI NG SEQUENCE : CALL CSIC 


PURPOS E ; Enter initial conditions on control system state 
variables into system state vector. 


COMM ONS USED ; ADS TAT, CSTAT, JCNTRL, IMAIN1 


CALLED B Y; MAIN 


SUBROUTINES CALLE D ; NONE 


GLOSSARY OF VARIABLES : See common block descriptions. 


DEBANG 


SUBROUTINE NAME: 


LAN GUAGE : FORTRAN IV 


CALLING SEQUENCE : 


CALL debang (yaw, roll, pitch) 


PURPOSE : Calculates the attitude of the spacecraft with re 

spect to the local vertical frame. The sequence 
used in going from the local vertical to the body 
frame is in the order yaw(3) , roll(l) pitch (2 ) . 

LABE LLED COMMONS USED: RP00L1, VI3CTRS 


CALLED BY : WHEEL S 


SU BROU TINES C ALLED : D SORT 


GLOS SARY OF VARIA BLES : OUTPUTS 

YAW ROTATION ABOUT THE THREE AXIS. 
ROLL ROTATION ABOUT TIIE ONE AXIS. 
PITCH ROTATION ABOUT THE TWO AXIS. 
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1 SUBROUTI NE NA ME ; DEB0UT 


LANGUAGE : FORTRAN XV 


CALLING SEQUE NCE: CALL DEB0UT 


PURPOSE : Provides the printed output for the control system 

l 

i using the subroutine SET. 

k 


LABE LLE D COMM ON S USED ; C STAT 


CALL E D BY : GPRINT 


SUBROUTINES . CAIiLED : SET 


GLOSS A RY O F VARI ABLES : See common bloclc description. 


■K ' 


SUBROUTINE! NAME ; DEREQ 


-t 


LANGUAGE : FORTRAN IV 


CALLING SEQUENCE ; CALL DEREQ 


PURPOSE ; Executive routine for derivative computation 


COMMONS USED ; See former DEREQ 1 


CALLED BY : DEREQ 1 


SUBROUTINES CALL ED : Sec former DEREQ 1 


GLOSSARY OF VARIABLES ; See former DEREQ. 


T CATI ONS : Eliminate argument;/, in calling sequence. DER] 

contains coding formerly in subroutine DEREQ 1. 


war*" 




. 
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SUB ROUT INE NAME i DJ3RKQ1 


LANG UAGE ; FORTRAN IV 


CALL ING SEQUENCE ; 


Call DEREQl (QDPEND, TIME, QDERIV ) 


PURPOSE ; Intermediary between subroutines ADMIMP and DEREQ 

DEREQ1 tabes state variables passed as an argument and 
places them in labelled common VARBLS . After calling 
DEREQ, this program also collects the calculated de- 
rivatives and returns them via the argument list. 


COMMONS USED ; IMAIN1 , RPOOLl. , VARBLS 

\ 

CALL E D BY : ADMIMP 


SUBR O UTINES CALLED : DEREQ 

GL OS SARY OF VAR IABLES : I NPUTS 

- ~ ~ ~ " TIME 

NUMEQS 

QDEPEND 


OUTPUTS 

TIM1 

DEPEND 

DERIV 


QDERIV 


Time in seconds from ADMIMP 
Number of state variables 
Values of state variables re- 
ceived from ADMIMP 


Time in seconds, passed to all 
other programs 

Values of state variables passed 
to all other programs 
Values of state vnrible deriva- 
tives received from computation 
programs 

Values of state variable deriva- 
tives returned to ADMIMP 


SUBROUTINE NAME : NUMCBE 


LANGUAGE : FORTRAN IV 


CALLING SEQUENCE : CALL NUMCSE 


PURPOSE: Determine number of control equations and defines mapp 

ing between system and control state vectors. 


COMMONS USED: ICNTRL , IMAINl, JCNTRL 


CALI, ED BY: NUM 


SUBROUTINES CALLED : NONE 


GLOSSARY OP VARIABLES: Sec common block descriptions 




SUBROUTINE NAME; PDTAPR 


LANGUAGE : FORTRAN XV 


CALLING SEQUENCE: Call PDTAPR 


PURPOSE: Begin reading ephemorio tape 


COMMONS USED: CNBODY, IMAINS , TJANl, VECTRS 


CALLED BY ; MAIN 


SUBROUTINES CALLED: DTAPRE , 1BCOM# 


GLOSSARY OF VARIABLES: See prior MAIN 


SUBROUTINE NAME : RWII0UT 


LANGUAGE : FORTRAN IV 

GALLING SEQUENCE: Call RWII0UT 

PURPOSE: Obtain trajectory data for external plot programs 

Writes to unit 50 when called. Called from ADMIMP 
after every time stop. 

COMMONS USED: CSTVAL, VARBLS 

CALLED BY: ADMIMP 


SUBROUTINES CALLED: II3COM# 


GLOSSARY OP VARIABLES: TTST Time in seconds 

DEPEND State variable values 


DERIV State variable derivatives 


SUBROUTINE NAME: VDIC 


LANGUAGE : FORTRAN IV 

CALLING SEQUENCE: Call VDIC 

PURPOSE: To compute initial conditions for viscous damper and 

insert into system state vector. 

COMMONS USED: ADSTAT , CONSTS, IMAIN1, RPOOL1, RVISCS 

CALLED BY; ‘MAIN 

SUBROUTINES CALLED; NONE 


GLOSSARY OF VARIABLES: 


See former MAIN program 


SUBROUTINE NAME: KXPN 


LANGUAGE : FORTRAN IV 


CALLING SEQUENCE: Call EXPK (T, TLA ST , TI VLAST, VI SSN, TCOR, 

OUT) 


PURPOSE: To generate noise channel-outputs for control system 

simulation . 

LABELLED COMMONS USED: ICSADM 


CALLED BY : WHEELS 


SUBROUTINES CALLED: GAUSS 


GLOSSARY OP VARIABLES: 


T 

TLAST 

TI 

VLAST (I) 
VI (I) 

SSN {A) 
TCOU(I) 
OUT (I) 


Current integrator time 
Time at last call to EXPN 


Time at last 
Noise channel 
Noise channel 
Noise channel 
Noise channel 
Noise channel 


successful integration time 
output at time TLAST 
output at time TI 
sigma 
lag 

output at time T 


3.4 niinnuPTioN or m;w ako modtftihv t,auei.i,kd common 

The new labelled common blocks are described in the follow- 
ing section. The only modification of labelled common in the 
0 OMMOK/Ti i id U5T/ • Tliis common wan changed to reflect the r.pecifi- 
cation of two pulses i>er spin period. 

COMMOM HI.OCK NAM1. ; 

COM-WNA'nRUPT/TVO^), TL0C ( 3 2), TTIM (4 , 2 ) , TPAR (4,2) , REF ( 2 ) 
The difi nitrons of thoso symbols are given in the input descrip- 


tion of Section 4.2 


£OMMOU 




NAMt 


COMMOi;/ADSTAT/l)IJU(X50) ( I>EP(150) 


USED . V. i SWMWUT:i } i 1 ;; l : MAIN, VDIC , CHIC 

PUR POSE ; The DER and DEP arrays \wre formerly local lo the MAIN 
Program. They have boon placed in common so that ini- 
tial condition can be loaded by VDIC, CSIC and future 
programs of the same general typo. 


VARIABLES : 

FORTR AN N AM? I TYPE 

DPI' (150) R*8 

DEP (150) R*8 


UNPl’S, DESCRIPTION 

variable State variables 


variable 


State variable 
derivative 


COfriMON/CSENDS/CSUP (20) , CSDN (20) , GNIC (10) 




* 

COMM ON ] Sll.-V IK NAME : 


USED IN SUBROUTINES: CSIC, KD478 


PURPOSE: To transfer integration bounds and noise model ini- 

tial conditions. 


VARIABLES : 

FORTRAN NAME TYPE 
CSUP (20) R*8 


CSDN (20) R*8 


GNIC (10) R*8 


UNITS DESCRIPTION 

varies Upper integration bound 
for control system state 
vector 

varies Lower integration bound 
for control system state 
vector 

varies Noise channel initial con- 
ditions 


* & 


COMMON RT r >GK NAME » COMMON/CN 0 ISE/VNS 2(10) , VNSI (10) , VNSN (10), 


Tl , T3 


USED IN SUBROUTINE S : CSIC, WHEEL S 


PURPOSE: Tran a 

for noiso 

channel 

data. 

VARIABLES : 




FORTRAN 

TYPE 

UNITS 

DESCRIPTION 

VNS2 (10) 

R*8 

varies 

Noise channel outputs for 
last call to EXPN . 

VNSI (10) 

R*8 

varies 

Noise channel outputs for 
last successful time step 
for integrator. 

VNSN (10) 

R*8 

varies 

Noise channel outputs for 
current call to EXPN. 

T1 

R*8 ' 

sec . 

Time of last successful 
integration 

T2 

R*8 

sec. 

Time of last call. 


** iifi'iiTr i 1 


COMMON BLOCK NAME: COMMON/C STAT/X (20), XDOT(20), CPARM(43) 


OR 


COMMON/CSTAT/SVCS (20) , SVDCS (20) , CPARM(43) 


U SED m S UBR OUTINES : RD47B, WHEELS’, CSIC, WRTPLT, DEB0UT 

/’ 


PURPOSE : Carries control system state vector, its derivative 

and control system parameters. 


VARIABLES ; 


FORTRAN NAME 

TYPE 

UNITS 

DESCRIPTION 

X 

R*S 

see 

Control system state 



inputs 

descriptions 

vector . 

XDOT 

R*8 


Derivative of control 
system state vector. 

CP ARM 

R*8 


Control system parameter 


,b. 




COMMON BLO C K NAME : COMMON/ICNTRL/KNT RL (10) 


USED IN S UBR OUTINES ; MAIN, RD470, WHEELS, NUM, NUMSCE, WRTPLT 


PURPO SE; Specifics control system options. 


VARIADJ.ES : 


FORTRAN NAME TYPE UNITS D ESCRIPTION 


KNTRL 


1*4 


See input descriptions 


COMMON iiI].p£K. „N AMJj; : COMMON/ICSADM/LSAVE, I RAND, NCI] AN 


USED IN S UBR OUTINES, i ADMXMP , CSIC, EXPN, WHEELS 


PUR POSE ; To transfer noise generator control words 


VAR TABLES : 


FORTRAN NAME 

TYPE 

UNITS 

DESCRIPTION 

LSAVE 

1*4 

N.A. 

Integrator logic integer 

IRAND 

1*4 

N.A. 

Random number .integer 


NCUAN 


1*4 


N.A 


Number of noise channels 


COMMON BLOCK NAME: COMMON/JCNTRL/NCNtRL t MCNTRL, MAPCNT(20) 


us E_D in . snimoim him ; csic, wheels, numcse 


i’tJIU w; : 


VAHXA BXiES ; 


FORTRAN NAME 

TYPE 

UNIT 

RECEIPT :r ON 

NCNTRL 

1*4 


Number of control equations 

MCNTRL 

1*4 


Index of state variable 
immediately proceeding 
1ST control state variable 

MAPCNT (20) 

1*4 


Defines mapping between 
system state vector and 
control state vector 


common/ sostat/ssvcs ( 20 ) 


CO MMON RIiO CK NAME : 


USED IN SUBRO UTINES ; CSIC 


PURPOSE ; To save control system state vector initial conditions 
for stacked cases. 


VARIA BLES : 

FORTRAN NAME 
SSVCS (20) 


TYPE UNITS 

R*8 varies 


DESCRIPTION 

Initial conditions for 
control system state 
vector 


COM MON B LOCK NAME : COMMON /VERS/VNO 


USED IN SUBROUTINES; 


MAIN, EVERS, (BLOCK DATA) 


PITR POBE : Identifies current program version used to determine 

modification level. 


VARIABLES : 

FORT R AN NAME TYPE UNITS DESCRIPTION 

R*8 


VNO 


8 character version 
identification 


4.0 DESCRIPTION OP NEW AND MODIP.TKD INPUT/OUTPUT 

The now and modified input symbols and output capabilities 
for the PSD program aro detailed in the following section. 

4.1 DYNAMICS EXPLOROR-R SIMULATION INPUT/OUTPUT 

The input symbols for the DE-B control system simulation 
are described below: 


INPUT PRESET 

SYMBOL TYPE VALUE DESCRIPTION OF USE 


KNTRL(IO) 1*4 0 Vector of control integers 

for DE-B control system 
simulation 


KNTRL 

ADDRESS DESCRIPTION 


( 1 ) 


( 2 ) 


(3) 

(4)-(8) 

(9) 


( 10 ) 


KNTRL (1)=0 No control system 
KNTRL (1)-1 Control system with 
second order sensors 


KNTRL ( 1 ) ~2 Control system with 
fourth order sensors 


KNTRL (2) =0 No nutation damper 
KNTRL (2) =1 Primary damper circuit 
KNTRL (2) "2 Offset pointing dm -or 
circuit 


KNTKL ( 3 ) — 0 No 

KNTRL ( 3 ) — 1 First order tachon., uer 


NOT USED 

KNTRL (9) -1 I is starting in- 
teger for noise gen- 
erator. I mu tit be 
odd and shoxxld have 
6 or 7 digits 

KNTKL ( 10 ) -0 No noise channels 

KNTRL (10) -3 Noise generated for 
sensors & bia.l volt- 
age. Use only 0 or 3 


¥ 


INPUT 

SYMBOL 

TYPE 

PRK53ET 

VALUE 

DESCRIPTION 

CPAPM (43) 

R*8 

(1-30) 0.000 
(31-40) 1.000 
(41-43) 0.000 

Control system parameters 
for DK-13 control system 
simulation 

CP ARM 
ADDREHP 

MATII 

SYMBOL 

UNITS 

DESCRIPTION 

(1) 

r s 

sec . 

Sensor time constant 

(2) 

r x 

see. 

Lead term in pitch compon^ 
sation 

(3) 

r 2 

sec. 

Lag term in pitch compen- 
sation 

(4) 

r F 

sec. 

Tachometer time constant 

(5) 

K s 

volts/rad. 

sensor gain 

(6) 

K c 

volts/volt 

Pitch amplifier gain 

(7) 

K a 

volts/volt 

Power amplifier gain 

(8) 

Kf 

volts/(rad/ see) 

Tachometer gain 

(9) 

K b 

volts/(rad/sec) 

Motor back EMF constant 

(10) 

K t 

ft. lbs/volt 

Morot torque constant 

(11) 



Not used for input 

(12) 

v lim 

volts 

VoH-.ayo limit in compen- 
satin’! amplifier 

(13) 

K /l 

vo3 ts/volt 

Ga. in primary damper 
circuit 

(14) 

% 

see . 

Time constant in primary 
dumper circuit 

(15) 

v b 

volts 

“Bias boltage 

(10) 

V'co 

ft-lbs 

Coulomb friction torque 

(17) 

-0- min 

* d/sec 

Tost relative wheel speed 



to avoid coulomb friction 
touquo discontinuity at 
nor ci speed 


CPARM 

ADDRESS 

MATII 

SYMBOL 

UNITS 

DESCRIPTION 

(18) 

So 

volts/volts 

Fourth order sensor 

(19) 

fy 2 

sec 

Gain in ofiVot pointing 
damper circuit 

(20) 

nrp 

N.D. 

Time constant in offset 
pointing damper circuit 

(21) 


N.D. 

Sign of damper circuit 
output set 10 1.0 or -1.0 

(22) 


volts 

Roll sensor* output limit 

(23) 


ft. -lbs. 

Motor torque output upper 
limit 

(24) 


ft. -lbs - 

Motor torque output lower 
limit 

(25-30) 



Not used 

(31-33) 



Noise model SIGMA for pitch 
roll and voltage bias re- 
spectively 

(34-35) 



Not used 

(36-38) 



Noise model LAG for pitch, 
roll and voltage bias re- 
spectively 

(39-40) 



Not used 

. (41) 


rad. 

Pitch sensor bias 

(42) 


rad. 

Roll sensor bias 


(43) 


Not tised 


INPUT 

SYMBOL 

TYPE 

PRESET 

VALUE 

SVCS(20) 

R*0 

0.000 

SVCS 

ADDRESS 


UNITS 

(1) 


volts 

(2-4) 


volts 

(5) 



(C) 


volts 

(7-9) 


volts 

(10) 



(11) 


volts 

(12) 



(13) 


volts 

(14) 



(15) 


rad/sec 

(1A-18) 



(19) 


volts 

(20) 


volts 


description 

Initial conditions for con 
trol system state vector 

DESCRIPTION 
Pitch sensor output 
Pitch sensor dynamics 
Not used 

Roll sensor output 
Roll sensor dynamics 
Not used 

State variable for pitch 
compensation amplifier 

Not used 

Tachometer output 

Not used 

Who el speed 

Not used 

Nu t a t ion damper 

Nutation damper 


INPUT 

SYMEOL 

GNIC(IO) 


CSUP (20) 


CSDN (20) 


PRESET 
TYPE VALUE 


DESCRIPTION 


R*B 0.0D0 


R*8 1.0D-2 


R*G 1.0D-4 


In5li.il conditions for noise 
model channels. 

GlITC(l) Pitch channel 
Gi;iC(2) Roll channel 
GK1C (3) Rian vollayo channel 
GNIC (4-1(1) Not unci 

Upper bound on difference 
between predicted and eorr- 
ecled control system state 
vector. Location in CSUP 
corresponds to the location 
of the variable in the 
state vector initial con- 
dition array SVCS 


Lower bound on difference 
between predicted and corr- 
ected control system state 
vector. Location in CSDN 
corresponds to the location 
of the variable in the state 
vector initial condition 
array SVCS 


In addition to the new input described above, the follow- 


ing must bn given s 
IWIIIUVO 1 


Forces call to momentum wheal subroutine 
(Wiim.f.) 


XMOM i b* ( 2 ) Inertia of momentum whe«*l (nlug/fi* ) 

It in recommended that integration control bounds bo given 
for critical control variables to prevent integration time stop 
from exceeding the time constants of the closed-loop system. For 
example, if fastest component in the control system has a time 
constant of .1 sec, it is unreasonable to expect accurate sim- 
ulation results with a larger time stop. Setting the inte- 
gration bounds to some reasonable (small) fraction of the nom- 
inal value (eg 10 to 10 ) will ensure that the integration 
errors are of the same order (eg 2 to 4 place accuracy) . The 
time stop corresponding to this level of accuracy will be det- 
ermined internally. 


The output for the DE-D control system simulation includes 
both printed data find plots. The printed output is as follows: 


OUTPUT 

HEADING UNITS 

PTCII Out 
ROLIi Out 
COMP Out 

TACII Out 

WHL SPD rad/soc 

NUTD Out 

The output available for plotting 

the control system. The location 

iables giro as follows; 

KPLOTS 

ADDRESS 

( 216 ) 

( 217 - 2 ] 9 ) 

( 220 ) 

( 221 ) 

( 222 - 224 ) 

( 225 ) 

( 226 ) 

( 227 ) 

( 22B) 


DESCRIPTION 

Pitch channel sensor output 

Roll channel sensor output 

Output of compensation 
amplifier 

Out|5Ut of tachometer 

Momentum wheel speed 

Nutation damper phase shift 
circuit output 

is the entire state vector for 
and definitions of these var- 

DESCRIPTION 
Pitch sensor output 
Pitch sensor dynamics 
Not used 

Roll sensor output 
Roll sensor dynamics 
Not used 

Output of pitch compensa- 
tion amplifier 

Not used 

Tachometer output 


ADDRESS 


UNITS 


DESCRIPTION 


(229) 

Nob used 


(230) 

Momon bum 

wheel speed 

(231-233) 

Not used 


(234) 

Nutation 

dynamics 

damper phase shift 

(235) 

Nutation 

output 

damper phase shift 


4.2 MODIFIED THRUSTER OPTION INPUT 


THRUST LO A DING OPTI ON 
FORTRAN MATH 

SYMBOL SYMBOL DESCRIPTION UNITS 


IPULSE N.A. 


ISPLSE 


N.A. 


ISPNP N.A. 


IPLPRP N.A. 


Control word to activ to N.D. 

IPULSE-0 No thrusting 
IPULSE 1 - 1 Apply thrust once 
IPULSE >1 apply thrust IPULSE 
times (only if ISPLSE 
* 0 .) 


Control word to activate sun N.D. 
crossing time to start the 
thrusting 


ISPLSE-0 Sun crossing not used 
I.SPLSE-1 Sun crossing used 

Control word to print out the 
orbit update message (only if 
ISPIiSE-'.l , IPULSE i) i.e., if 
ISPNP-5, the orbit update mess- 
age will bo printed at every 
?.*th pusla 

Control word for number of N.D. 
thrust pulses per spin record 

IPLPRP-* 1 One pulse 
IPLPRP -2 Two pulses 

Only one or two pulses arc 
allowed, (Preset : 1) 


FORTRAN 

SYMBOL 

TVKCTR (3,2) 

TLOCAT (3,2) 
Tl’IMES (4,2) 


TPARAM ( 4 , 2) 


MATH 

SYMBOL DESCRIPTION UNITS 


(By) 



t. 


t 


4 


A 

B 


C 

D 


Unit vector defining the, N.D. 

direct! on of i he force 
applied to the. body due 
to thrusting . This vec- 
tor is defined in the 
body frame . (preset^ 

0 . 0 , 0 . 0 , 1 . 0 ) 

Location in the body feet 

frame of the point app- 
lication of the force 
due to thrusting. 

(proset-0.0) 


Times to define thrust var- 
iation measured from the 
problem starting time 

TT1MES (1,1) Start of pulse 

TTXMES (2,1) End of exponen- 
tial rise 


TTIMRS (3,1) End of linear 
thrust 

TTXMES (4,1) End of pulse 
1=1 or 2 


Parameters to define thrust 
variation 


sec 


TPARAM (1,1) Coefficient dur- lb 
ing exponential rise 

TPARAM (2,1) Exponential decay sec 

constant durint dxponential 

rise 


TPARAM (3,1) Coefficient for lb/sec 

linear slope 

TP ARAM ( <1,I) Exponential decay sec ^ 

constant during exponential 

decay 


FORTRAN 

SYMBOL 

MATH 

SYMBOL 

DESCRIPTION 

UNITS 

RNFANG ('.».) 

A, 

U 

Angular dolay from tho Y^ 
ax.ir. erofu?.ing tho sun lane 
to tin 1 pul. to 

dog 

4.1 MODIFIED 1 

*F i' INPUT AK» 

ADD IT TONAL PLOT CAPABILITY 


ldi nt, Fpur.i < 'iy J f 'ji 

ynpfqrm.iFFT) 

..Anglyri r 


FORTRAN 

SYMBOL 

MATH 

SYMBOL 

DESCRIPT1 ON 

UNITS 

XCSD 

N.A. 

Control word to activate the? 
FFT analyst's .subroutine. 

N.I). 


ICSlt-0 No FFT anal ysir. 

ICSD-N FFT analysis activated 
(pr 0:30 h ; 0) 

Tho integer N requests that N 
frequencies bo extracted from 
the? .selected I< PLOTS data .sols 

NOTE; Tho value of’ N should not oxeeod tho number of froquon- 
cior, that can reasonably bo expected to exist in tho 
data. Tho range of permisfublo valuer. for N art' from 1 
to 10. 

1 ^ N 10 




* atm**.--.**.****.- 




« 




KPEOTS 

ARRAY 


address 

i 

DES0R1 1'T.rON 

FORTRAN 

SYMHOTi 

UNITS 

NOTE 

:no 

External moment 
about 1 body axis 

MOMEN I’ 1 

ft -lbs 


211 

Ext ornal moment 
about 2 body axis; 

MOMENT 2 

ft- lbs 


PI o 

u Jj M 

External moment 
about 3 body axis 

MOMENT 3 

ft- lbs 


213 

Component of ang- 
ular momentum on 
1 body ax if? 

IIRODY 1 

s.l ug--£t "'/see 

incAiiO i 

214 

Componant of ang- 
ular momenlun on 
2 body axis 

I1H0DY 2 

p 

slug-f t" dsoc 

XllCATiC 1 

215 

Componant of ang- 
ular momentum of 
3 body axis 

111 >04 it 3 

slug/f ti’/sec 

.1 IK' Alp 1 

21b 

t 

Coni rol system 
state vector 
(See Section 4.1) 

N.A. 

varies 

IWllKEb 1 
and 

RNTRI.(l) 


1. or 


V 


